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Abstract Recently Pop (Solar Phys. 276, 351, 2012) identified a Laplace (or 
double exponential) distribution in the number of days with a given absolute 
value in the change over a day, in sunspot number, for days on which the sunspot 
number does change. We show this phenomenological rule has a physical origin 
attributable to sunspot formation, evolution, and decay, rather than being due to 
the changes in sunspot number caused by groups rotating onto and off the visible 
disc. We also demonstrate a simp l e met hod to simulate daily sunspot numbers 



over a solar cycle using the IPopI (|2012|) result, together with a model for the 
cycle variation in the mean sunspot number. The procedure is applied to three 
recent solar cycles. We check that the simulated sunspot numbers reproduce the 
observed distribution of daily changes over those cycles. 

O J Keywords: Solar Cycle, Models; Sunspots, Statistics 

1. Introduction 

Sunspots have been studied scientifically since the invention of the telescope, 
and reliable daily sunspot number records are available from the early 180ufl 
The accepted quantity to characterise solar activity is the international sunspot 
number 

s = k(10g + n), (1) 
where n is the number of indi vidual sunspots, q is the num ber of sunspot groups, 



and A: is a correction factor (jBruzek and Durrantl . 119771) . The sunspot number 



changes in a secular or long-term fashion with the semi- regular 11-year sunspot 
cycle, driven by an internal magnetic dynamo which generates the magnetic 
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field (jTobiasl . 120021) . The secular change in sunspot number exhibits apparent 
randomness, as evidenced by the extensive literature describing the smoothed 
daily sunspot number (typically a monthly average) as a stochastic, or chaotic 
time series. In particul ar, the peak and timing of the different solar cycles shows 
considerable variation (|PetrovavLl2010h . 

Sunspot numbers also vary on shorter time scales, in particular daily, as a 
result of the complicated local processes associated with sunspot formation, 
evolution, and decay. The short-timescale variation pr oduces large excursions 
in su nspot number, up to 100 per day in extreme cases ( Noble and Wheatland! . 
20111 ). The large day to day variations are caused by the rapid evolution of mag- 
netic structures, and the sudden appearance/disappearance of large active re- 
gions. These rapid developme nts can have important consequences for the space 

weat her experienced on Earth (jCommittee on the Societal and Economic Impacts of Severe Space Weather I 
20081) . 

Recently, it was demonstrated that the change As in the daily sunspot num- 
ber (for days on which the number does change) follows a Laplace, or double 
exponentia l distribution (jPopl . 120121 ). Including the case of days with no changes, 
Pod ( 2012f ) modelled the distribution /(As) of the change in sunspot number 
with the functional form 



/(As) = AI(As < 0) exp(As/B) + AI(As > 0) cxp(-As/B) + CS(As), (2) 

where A, B, and C are constants, and where I(x) is the indicator function defined 
by I(x) = 1 for x true, and I(x) = for x false. The parameter C determines the 
fraction of zero changes, and B is the mean absolute change for days on which the 
number does change. Normalisation of Equation @ requires A = (1 — C)/2B. 

Figure [T] shows the observed distribution for the daily sunspot number using 
th e NGDC da ta for 1850-2011, and illustrates the exponential form identified 
by IPqpI ( 20121 ). The top panel is a histogram of daily changes As. The bottom 
panel shows the cumulative number of changes greater than Asj, for positive 
changes: 

N{As > Asi) = K As > As i > 0), (3) 

i 

and the cumulative number less than Asi for negative changes: 



iV(As < As,) = ^I(As < As, < 0). 



(4) 



Both panels use a logarithmic scaling on the vertical axis. The red dots in the 
panels show the data, and the blue curves show the model distribution defined by 
Equation (ffl). with param eters B and C estimated from the data using maximum 
likelihood (|Ehasonl . 1 19931 ). The adherence to the exponential form in Figure [T] is 
strik i ng. T here are a large number of days with no change in sunspot number, and 
Pol (|2012 l) refers to As = a s a "sp ecial state" . The distribution is remarkably 
symmetric about As = 0. Pop! ( 2012 ) also investigated the solar cycle dependence 
of the observed distribution, and found that the exponential rule in absolute 
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changes is most closely adhered to over whole cycles, with the greatest departure 
near minima of cycles. 

It is surprising that the Laplace distribution in the change in daily sunspot 
number was not identified and discusse d in the literature earlier. The behaviour 
was previously noted in smoothed data (jLepreti. Kossobokov. and Carbond . l2009h .| 



and an approximate exponential dependence in the distrib ution of overall sunspot! 



numb ers, related to Equatio n ((2]), was also commented on Noble and Wheatland! 



( 201l[ ). However, Pop ( 2012 ) showed that the adherence of the observed changes 



in daily sunspot number to the Laplace distribution defined by Equation @ is 
much stricter than the approximate exponential distribution of overall sunspot 
number. The Laplace distribution (Equation @) represents a newly identified 
phcnomcnological rule describing the way in which the daily sunspot number 
varies, which should have application for modelling and prediction. The origin 
of the distribution is not obvious. Changes in sunspot number occur daily due 
to sunspot group formation, evolution and decay, the appearance/disappearance 
of spots and spot splitting (referred to here as spot evolution), and also due to 
sunspot groups and individual spots rotating onto and off the visible disc. If the 
law is due to spot and group evolution, then the phenomenological rule must 
have a physical origin. 

In this paper we demonstrate that the Laplace distribution of changes in 
sunspot number is caused by sunspots forming, evolving, and decaying, and 
is not a result of rotation on and off the disc. We then also present a simple 
Monte Carlo method, based on Equation © and some additional assumptions, 
for simulating sunspot numbers over solar cycles. 



2. Origin of the Laplace Distribution 

To investigate the origin of the observed distribution in changes in sunspot 
number we use reports of sunspot groups on the Sun, for 1981-2011, compiled 
by US National Oceanic and Atmospheric Administration (USAF/NOAA)OThe 
daily change in sunspot number 

A.s = fc(10A.g + An) (5) 

can be decomposed into changes due to rotation of regions and spots onto and 
off the disc As r , and changes due to group and sunspot evolution As e : 

As = As r + As c . (6) 

Similarly the terms Ag and An on the right hand side of Equation (|5|) can be 
decomposed in this way. To approximate the change due to rotation As r we 
assume that active regions first appearing on the disc within 180/14 w 13° of 
the eastern limb arrive due to rotation within a day, and regions last observed 
on the disc within 13° of the western limb disappear due to rotation within a 



2 Data are available at http://ngdc.noaa.gov/stp/solar/sunspotregionsdata.html 
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day. The factor of 14 days is the approximate time to rotate across the disc 
( Snoderass and Ulrich . 1990t) . 

Figure [2] shows the result of the analysis of the data. The figure presents 
the cumulative distribution of the total change in sunspot number As using the 
complete data set (black), the change due to rotation of spots and groups onto 
and off the disc As r (blue), and the changes due to the evolution of spots and 
groups As c (red) . The distributions of the total change As and change due to 
evolution As e are very similar, and both clearly show the Laplace distribution 
of Equation ^ . The distribution of As r exhibits significant falls in number at 
As r = ±10fc, ±20fc, which may be attributed to entire sunspot groups rotating 
onto and off the disc (groups are weighted with a factor of 10 in the definition 
of the international sunspot number - see Equation (|T))). 

An important feature of Figurc[5]is that the number of changes due to rotation 
onto and off the visible disc is at least an order of magnitude smaller than 
the number of changes due to evolution. There are too few changes associated 
with active regions rotating onto and off the visible disc for this to influence 
the overall distribution of changes in sunspot number. It is very unlikely that 
the observed double exponential distribution of changes, for days on which the 
sunspot number does change, is influenced by the small number of active regions 
which rotate onto or off the disc in a day. 

It is possible that the observed Laplace distribution arises from an averaging 
of locally non-exponential distributions across the disk. To test this we also 
calculate the cumulative fraction of changes occurring in particular 26° strips 
on the disc. In other words, we identify the active regions in the strips on a 
given day, and calculate the sunspot number due to these regions only. On the 
next day we identify the regions in the strips, and again count the sunspot 
number due to these regions. The difference between these two counts is the 
change in sunspot number in a day occurring in the particular strips. This 
count includes both changes due to rotations into and out of the strips, and 
changes due to evolution in the strips. Figure |3] shows the results. Changes 
occurring in the region (— 13°, 0°) and (0,13°) are in blue, changes occurring 
in the region (—35°, —22°) and (22,35°) are in black, changes occurring in the 
region (—57°,— 44°) and (44,57°) are in green, and changes occurring at the 
limbs (i.e. (-90°, -77°) and (77,90°)) are in red. Figure shows that changes 
associated with regions in a restricted range of longitudes reproduce the observed 
Laplace distribution. The exception is the distribution of changes at each limb, 
which shows large falls in probability due to changes in the number of groups 
in the strips, and a significantly smaller proportion of large changes in sunspot 
number. Presumably these differences arise due to the difficulty of accurately 
observing sunspot regions at the limb, and in particular, resolving the number 
of individual spots in a group. This may explain why this particular distribution 
is dominated by changes in the number of groups. 

Figure [3] shows that the observed Laplace distribution in total changes does 
not arise from averaging over non-exponential local distributions in particular 
strips. The conclusion is that the Laplace distribution form is representative of 
changes in daily sunspot number due to the local evolution of sunspot groups 
on the Sun. It reflects the physical processes occurring. 
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3. Daily Change in Sunspot Number 

3.1. A Conditional Distribution for Daily Change in Sunspot Number 



Pod ( 2012 ) investigated the daily change in sunspot number As using data over 



complete cycles, which involve days with a range of different initial values s 
of the sunspot number at each change. The phenomcnological rule represented 
by Equation (0) was found to hold to a very good approximation for the range 
10 < As < 60 over the last 14 cycles, although departures from the rule for small 
changes in sunspot number were noted. This departure may be attributed to 
the discrete nature of sunspot number, which means that the minimum sunspot 
number larger than zero is 11 A: (where k is typically less than unity). This causes 
a discrete jump in the tabulated daily values of the international sunspot number 
from zero to seven (and implies that the average value of k used by observers is 
k = 0.64). Departures for large changes in sunspot number were also noted. This 
may be attributed to the finite size of the sunspot number over any cycle. Large 
negative changes in sunspot number are unlikely because the sunspot number is 
unlikely to have a sufficiently large value at any given time over a cycle to allow 
that change. 

The distribution of a change As on a given day is dependent on the value of 
sunspot number on that day. To model this we introduce transition probabilities 

p(s -)• s'\s) =p(s',s) (7) 

for changes from an initial sunspot number s to a final sunspot number s' on 
a day, given that the sunspot number is initially s. The function on the right 
hand side is a conditional distribution, and in Section [3.2l we consider a suitable 
functional form for this distribution. In Section 13.31 we relate the chosen form 



of the conditional distribution p(s' , s) and I Pod 's exponential distribution /(As) 
of changes over complete cycles, and demonstrate how to estimate parameters 
from the data, for the proposed conditional distribution. 



3.2. The Form of the Conditional Distribution 



To gain insight into a suitable function form for the conditional distribution 
in Equation ([7]) we re-examine the data. Figure [4] is a two-dimensional (2-D) 
histogram A^(As, s) of the number of days for which the sunspot number simul- 
taneously has the value s (enumerated along the vertical axis) , and increases by 
As to the next day (horizontal axis), for the NGDC data for 1815-2010. The bins 
are chosen to be of size two in As and s, and the figure shows the normalised 
histogram 

p itj = N(A Sll Sj )/ ^ A(As,, Sj ), (8) 

i 

so that the greyscale density along each row shows the relative probability of a 
given change As^, for the given initial sunspot number Sj. A nonlinear scaling is 
applied to the greyscale density, to better show the bins with small numbers 
of days. Figure 0] illustrates the influence of the lower boundary As = — s 
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required by the non-negativity of sunspot number. The distribution is relatively 
symmetric about As = for any given initial sunspot number s, except for an 
excess of days at the As = — s boundary (which correspond to changes leading 
to a zero sunspot number), and an excess of days with no change (i.e. As = 0). 

Based on Figure 21 a suitable approximate model form for Equation ([7]) is a 
simple exponential distribution symmetric about As = for each value of s, with 
the same coefficient in the exponent for both negative and positive changes. The 
(asymmetric) non-negativity of sunspot number may be imposed by requiring 
that transitions producing a negative final sunspot number (s' < 0) lead to zero 
sunspot number (s' = 0) instead. This may be written as: 

p(s',s) = Dexp[-(s - s')/E]l(s' < s) + L>cxp[(s - s')/E]l{s' > s) 

+GS(s') + HS(s' -s), (9) 

where D, E, F, G, and H are constants, and where the two terms involving delta 
functions describe the excess of days with changes leading to zero final sunspot 
number, and the excess of days on which the sunspot number does not change, 
respectively. The data show that for the observed changes in daily sunspot num- 
ber for 1850-2011 the fraction of negative and positive jumps are approximately 
equal (42% and 41% respectively), so we assume the same symmetry holds for 
each value of s in Equation 



o 



p(s',s)ds' = J p(s',s)ds'. (10) 
Normalising over all final sunspot numbers s', i.e. requiring 

p(s',s)ds' = l (11) 



leads to 



and 



G = 1(1 - H)e-' E . (13) 

The model conditional distribution (Equation ©) then has two parameters, E 
and H . The parameter E determines the size of changes on days when there is 
a change (a typical value, based on the data, is E rj 10). The parameter H is 
the probability of no change in daily sunspot number. 



3.3. Relating the Conditional Distribution and the IPopI (|2012l ) Distribution, and 
Parameter Estimation 

The overall distribution of changes /(As) may be calculated from the transi- 
tional distribution (Equation (0)) by integrating over all starting values s, i.e. 
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calculating 

/>oo 

f(As)=L p(s + As >S )g(s)ds, (14) 
Jo 

where L is a constant imposing normalisation over changes in sunspot number: 

/oo 
/(As)dAs = 1, (15) 
-oo 

and g(s) is the probability of an initial sunspot number s in the observations. 
A suitable choice to approximately describe the distribution of sunspot numbe r 
g(s) over a complete solar cycle is an exponential ( Noble and Wheatland , 2011 ^1 



g(s) = i e -/\ (16) 

where the mean value of daily sunspot number (based on observations for 1850- 
2011) is A = 55. Calculating the integral in Equation (fT4"f gives 

Lf(As) = ^£-c As ' E [1 + f e As / A ] I(As < 0) + ^- C - As / E I{As > 0) 

+H5(As) (17) 

and the normalisation by Equation (|15[) implies 

E(l-H) , , 

For the case of negative changes (As < 0), Equation (JTTJ) contains a term 
A -1 Ee As / x which makes the distribution /(As) asymmetric about As = 0, and 
which is produced by changes leading to zero sunspot number. The characteristic 
size of this term is E/X « 0.18, suggesting that the asymmetry in the distribution 
in As produced by the lower boundary As = — s is not a strong effect. In the 
absence of the extra term the normalisation constant is L = 1, and based on 
the data we find L w 1.05. Neglecting the extra term and setting L = 1, the 
distribution of changes implied by the conditional distribution (Equation ^) is 

/(As) = ±=£e As / E I(As < 0) + ^fc- As / E I(As > 0) + HS(As) (19) 



which is the same functional form as Equation ©, the lPod (|2012h distribution 
of changes. 

The correspondence between the conditional distribution (Equation @) and 
the overall distribution of changes (Equation (fTO"]) ) allows the parameters E and 
H of the conditional distribution to be estimated from the daily changes As over 
a cycle using maximum likelihood. Specifically, the exponential coefficient B (see 
Equation ©) estimated for the overall distribution may be taken as the estimate 



3 A noted in Section [2] the adherence of the overall sunspot number to an exponential distri- 
bution is only very approximate (by comparison with the observed distribution of changes in 
daily sunspot number, which follows the Laplace distribution quite strictly). 
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for coefficient E in the conditional distribution, and the fraction of days C with 
no change in sunspot number (see Equation @) in the overall distribution may 
be taken as the estimate of the corresponding parameter H in the conditional 
distribution. 



4. Simulating Sunspot Numbers 

In this section we apply the conditional distribution (Equation (jH])) in a Monte 
Carlo simulation of daily sunspot numbers Sj = s(tj), where ti refers to a day. 
The daily random variation in sunspot number is described by the stochastic 
differential equation (stochastic DE) 

ds 



dt 

where ./V is the number of days and the daily change Asi is generated by sampling 
from the distribution 

p(s' i , Si) = p(si + As^ Si) (21) 

with p(s\ s) given by Equation ©. 

To solve the stochastic DE (Equation (|2TJ)) ) we need to sample from the condi- 
tional distribution p(s', s) given by Equation ©, with our maximum likelihood 
estimate of the parameters E and H. This is achieved as follows. For a given 
initial sunspot number s and the estimates of E and H we generate a final 
sunspot number s' = s + As by generating an exponential random variable As 
with parameter E i.e. a random deviate As distributed according to ~ e~ As / E . 
A random variable u which is uniformly distribution on (0, 1) is also generated, 
and then the final change As is calculated according to the rule: 

• if u < 0.5(1 - H), then As = -As: 

• if u > 0.5(1 +H), then As = As; 

• otherwise As = 0. 

This procedure assigns no change in sunspot number (i.e. As = 0) with probabil- 
ity H, and the remaining changes are exponentially distributed over positive and 
negative As with equal total probability. The mean absolute size of changes (on 
days when As ^ 0) is E. Finally, to prevent the final sunspot number s' = s + As 
from being negative, if s + As < 0, then we take As = — s. 

This procedure assigns values correctly according to Equation Equa- 
tion ([2TJ]) is solved for a given initial sunspot number sq at time to by generating 
a sequence of transitions As^ (with i = 1,2,..., N) according to this recipe, and 
adding these successively to s . 

This model accounts for the stochastic variation in sunspot number according 
to the conditional distribution given by Equation © , but it does not account for 
the secular or long time-scale variation of the sunspot number over a solar cycle 
(the "shape" of the cycle). Recently iNoble and Wheatland! ( 2011 ) presented 
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a method for modelling the solar cycle variation in a general Fokker-Planck 
description of stochastic variation in sunspot number, and the same procedure 
is applied here. A term may be added to the stochastic DE (Equation ([2"0]) ) 
causing the fluctuating sunspot number to return to a prescribed time evolution 
9{t): 

d N 

JL = K {9(t)-s}+Y,^- (22) 

i=l 

The function 9(t) is referred to as the driver function, and factor k is the rate 
at which sunspot number s returns to the value specified by the driver function. 
The two terms on the right hand side of Equation ([22]) are deterministic and 
stochastic terms, respectively. Equation (|22[) is solved for a given initial sunspot 
number sq by adding daily stochastic transitions As^ in the same way as for 
Equation (|20p . In between the transitions the sunspot number is evolved accord- 
ing to Equation (|22p with just the deterministic term included. The solution to 
the differential equation with just the deterministic term is 

s*(t) = e -<*-**) Li + k f 9{t')e Kt ' dt'X , (23) 



where Sj = s(^) is the value of the sunspot number on the most recent day, and 
ti < t < 

To summarise, the procedure for simulating the sunspot number evolution for 
day i + 1, given the sunspot number Sj on day i, is to evaluate a deterministic 
value for the sunspot number s*(ii+i) using Equation (|23p . to generate a random 
change As^+i using Equation ((9|) with initial sunspot number s = ,s*(i; + i), and 
then the sunspot number on day i + 1 is Sj+i = s*(tj+i) + Asj+i. 

The driver function 9{t) in the deterministic term in Equation (|22[) represents 
the functional form of the solar cycle variation in sunspot number, i.e. the shape 
of a cycle. The driver function describes basic empirical features of a solar cycle, 
such as the time taken to reach maximum, the size of the maximum, and so on. 
Additionally, 9(t) may a ccount for detailed features of the shape of a cycles such 
as the Gnevyshev Gap (lGnevvshevl.ll967l). Here we use a function introduced 



by lHathawav. Wilson, and Reichmannl (jl994l hereafter HWR94) : 



exp [— [t — toy /b z \ — c 

where to is the start time for a cycle, and a, b, and c represent the cycle 
amplitude, period, and asym metry, respectively. Equation (j24|) was used by 
Noble and Wheatland! ( 2012 ) to investigate daily variation in sunspot number, 



and to forecast sunspot number and solar cycles. A statistical procedure for cal 
culating estimates of the parameters a, 6, c, and n from daily sunspot data, and 



the va lues of the estimates for cycles 11-23, was given in lNoble and Wheatland 



r)2012h . Here we re-use these parameter estimates describing the shapes of the 



cycle to simulate three recent solar cycles (21, 22, and 23). 
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Figure [5] shows the daily sunspot numbers over cycle 23, for the years 1996 
to 2008 (red points), and our simulation of sunspot numbers for this cycle based 
on Equation (|2"2"|) (blue points). The HWR94 driver function, Equation ({Ml) , 
enforces the secular variation in the solar cycle, with the parameter values 
a = 7.82 x IP" 8 , b = 15 1 4, c = 0.222, and K = 0.086 (taken from Table 
I in Noble and Wheatland (|2012h ). We also use the estimates E = 8.78 and 



H = 0.149 for the conditional distribution (Equation ([9])), which are maximum 
likelihood values on the changes in daily sunspot data for cycle 23, as discussed 
in Section [3] For the HWR94 driver function it is not possible to solve Equation 
([23)1 analytically, and instead we integrate 

J = *[*(*)-»] (25) 



numerically using a fourth order Runge-Kutta scheme ( Press et a?J . ll992l ) 



Figure El shows the daily changes in sunspot number for cycle 23 (red points), 
and the corresponding changes in our simulation (blue points) . The upper panel 
of Figure [5] shows the distribution of changes and the lower panel shows the 
cumulative distribution in the same format as Figure [1] Figure [6] confirms that 
the simulation of daily sunspot number over cycle 23 based on the stochastic 
DE (Equation (|22p) and the conditional distribution (Equation ©), together 
with the HWR94 model for the shape of the solar cycle, generates a distribution 
of changes /(As) over the cycle that closely resembles the exponential form 
identified bvlPopI (|2012t) . 



Figure [7] presents the results of the simulation procedure for cycle 22 (years 
1986 to 1996) in the same format as Figure [6] The parameter estimates for the 
HWR94 driver function are a = 14.1 x 10~ 8 , b = 1368, c = 0.33, and k = 0.073, 
again taken from Table I in Noble and Wheatland! ( 2012 ). The estimates E = 



10.4 and H = 0.096 arc used for the parameters in the conditional distribution 
Equation (J^J) , based on maximum likelihood applied to the daily data for cycle 
22. 

Figure |S] presents the results of the simulation procedure for cycle 21 (years 
1976 to 1986), again in the same format as Figure [6l The parameter estimates 
for the HWR94 driver function arc a = 12.2 x 10~ 8 , b= 1414, c = 0.490, and n 



0.073, again taken from Table I in lNoble and Wheatland! (|2012h . The maximum 



likelihood estimates E = 10.6 and H = 0.098 are used for the parameters in the 
conditional distribution (Equation 

Figurcs[7]and[5]confirm that the simulation procedure succeeds in reproducing 
the phcnomcnological exponential rule for daily changes in sunspot number when 
applied to cycles 22 and 21, respectively. 



5. Discussion 

This paper establishes that the observed Laplace, or double exponential distri- 
bution of changes As in daily sunspot nu mber s (for days on which the sunspot 
number does change) recently identified by Porj ( 2012t ). is due to the evolution of 



observed sunspot groups (i.e. group formation, spot splitting, spot/group decay) 
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rather than being due to the artificial variation caused by groups rotating onto 
and off the visible disc. The implication is that the distribution has a physical 
basis. Sunspot emergence, evolution, and eventual decay produces daily changes 
in sunspot number which may be positive or negative, and changes of this kind 
in separate active regions may add or cancel. The sum of these daily changes, 
remarkably, produces a simple Laplace distribution, with a marked excess of 
days with no change in sunspot number. 

In this paper we show also how to simulate daily sunspot number via a Monte 
Carlo method, using a conditional distribution based on the exponential rule 
together with a model for the solar cycle variation in sunspot number. The 
conditional distribution p(s', s) introduced describes the probability of a change 
from a current sunspot number s to a value s' = s + As in one day, given 
the initial sunspot number, and ensures that s' > 0. The simulation procedure 
involves calculating a secular or deterministic change in sunspot number due 
to the underlying solar cycle, and then adding a random change in sunspot 
number according to the conditional distribution. The Monte Carlo method is 
demonstrated in application to three recent solar cycles (cycles 21, 22, and 23). 
The simulated sunspot numbers exhibit a distribution of changes /(As) over 
each c ycle that closely reproduces the Laplace distribution identified by 



(|2012h 



It is interesting to consider possible explanations for the observed Laplace dis- 
tribution. The origin of the surface changes in sunspot number described by the 
rule are changes in the structure of the su bphotospheric magnetic fi elds, which 
are not directly amenable to observation (jThomas and WeissL Il99lf). although 



local helioseismology is beginning to provide some insights ( Gizon and Birchl . 



20051 ). In the absence of detailed physical models for the surface changes pro 



vided by this field evolution, it may be possible to construct statistical models 
for daily changes in sunspot number based on simple statistical descriptions of 
spot and group evolution, for given numbers of spots and groups. For example, 
probabilities could be assigned to given spots or groups increasing or decreasing 
their number in a day. Such a description may be modelled by a type of birth- 
death process, which have been used in the natural and social sciences (e.g. 
see for example iGillespid . Il992h . It may also be possible to use other known 



statistical rules for the distribution and evolution of spot gr oups, for example 



spot gr c 

the log-normal distribution of spot areas (|Bogdan et al.lll988l) . and various rules 
for the de cay rate of sun spot area (in area per unit time) per sunspot within a 
group (see SolankiL 2003 ) . Given the complexity of the combinations implied by 



this consideration of possible modelling, we note again how remarkable it is that 
a simple double exponential form arises. We leave the development of a detailed 
model to a future investigation. 
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